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ABSTRACT 


Using the low-frequency cutoff of electromagnetic noise trapped 
in the magnetosphere at frequencies abo/e the local plasma frequency 
it is now possible to make very accurate, + 1$, electron density 
measurements in the low density region between the magnetopause and 
plasmapause. This technique for measuring the total plasma density 
has been used, together with measurements of the suprathermal proton 
intensities with the LEPEDEA instrumentation on the IMP-6 spacecraft, 
to determine the thermal proton densities in the region between the 
plasmapause and magnetopause . Although the thermal protons usually 
account far a significant fraction, ~ 50$, of the total proton 
density in this region, in some cases, particularly at the larger 
radial distances the density of the thermal protons sometimes drops 
to a very small fraction, < 5$, of the total density and nearly all 
of the plasma consists of suprathermal particles. 

Usually the high-energy tail of the thermal proton distribution 
is sufficiently intense to be detected at low energies within the 
LEPEDEA energy range. When the high-energy tail of the thermal 
distribution can be detected, it is possible to determine the tempera- 
ture of the thermal protons. The temperature in a typical case is 
found to be 80,000°K. This high temperature for the thermal protons 
which are presumably of ionospheric origin agrees with earlier 0G0-5 
measurements reported by Serbu and Maier. 
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I . INTRODUCTION 


The ratio of the cold thermal plasma density to the hot supra- 
thermal plasma density plays a crucial role in trie stability of 
magnet ospheric plasmas [Kennel and Petschek, 1965 ; Cornwall et al. , 
1970; Brice and Lucas, 1971]* Qualitative comparisons of the thermal 
and suprathermal plasma densities have been made by Russell and 
Thorne [1970], using observations reported by Frank [1967a] and 
Taylor et al. , [ 19b8] . These data show that the hot ring-current 
protons penetrate deep into the cold ionospheric plasma inside the 
plasmapause during a geomagnetic storm. The region of overlap between 
the cold ionospheric plasma and the hot ring-current plasma is 
unstable far ion-cyclotron waves which resonate with the ring-current 
protons [Cornwall et al. , 1970] . Similarly, the penetration of supra- 
thermal electrons into the plasmasphere is believed to cause plasma- 
spheric hiss, via the resonant interaction of these electrons with 
whistler -mode waves [Thorne et al. , 1973]* In both cases the presence 
of the cold plasma contributes significantly to the instability by 
lowering the energy at which resonance occurs, thereby increasing the 
number of particles which are in resonance with the wave. 

Although there have been numerous independent measurements 
of both the cold plasma [Taylor et ai. , 1968; Serbu and Maier, 1970; 
Chappell et al., 1971] and the hot suprathermal plasma [Frank, 1967b; 
Vasyliunas, 1968; DeForest and Mcllwain, 1971] there have been no 
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simultaneous measurements of both the thermal and suprathermal 
plasma densities in the low density region between the plasmapause 
and magnetopause with sufficient accuracy to be useful for quanti- 
tative investigations of the plasma instabilities which occur in 
this region. The purpose of this paper is to present a series of 
total electron density measurements using data from the University 
of Iowa plasma wave experiment on IMP-6, and simultaneous measure- 
ments of the suprathermal, 52 eV < E < 38*000 eV, proton densities 
with the LEPEDEA instrumentation on IMP-6. With these measurements 
we are able to determine the density and temperature the thermal 
protons and the relative contributions of the thermal and suprathermal 
components to the total proton densities in the region between the 
plasmapause and magnetopause. 
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II . INSTRUMENTATION 


The plasma wave and charged particle measurements reported 
here were obtained with instrumentation on the satellite IMP-6 
(Explorer 43). IMP-6 was launched on March 13, 1971, into an orbit 
with initial perigee and apogee geocentric radial distances of 
6620 km and 211,250 km, respectively, and an inclination of 28.6°. 

The plasma wave antennas on IMP-6 comprise two long dipole antennas 
with tip-to-tip lengths of 53*5 and 92. 5 meters for electric field 
measurements and three mutually orthogonal loop antennas for 
magnetic field measurements. Two 16-channel spectrum analyzers, 
one connected to an electric antenna and the other connected to a 
magnetic antenna, are used to make spectrum measurements over the 
frequency range from 20 Hz to 200 Hz. Two wideband receivers are 
also included to provide broadband coverage of the frequency ranges 
from 10 Hz to 1 kHz and from 650 Hz to 30 kHz, depending on the 
particular mode of operation selected. Further details of the plasma 
wave instrumentation on IMP-6 are given by Gurnett and Shaw [1973]. 

The Low Energy Proton and Electron Differential Energy 
Analyzer, or LEPEDEA, on IMP-6 provides measurements of the differ- 
ential energy spectra and angular distributions of proton and 
electron intensities over the energy range 10 eV < E <, 40,000 eV. 

These electro? tic analyzers are capable of simultaneous determinations 



6 


of proton and electron intensities within each of 16 energy band- 
passes. The axes of the fields -of -view of the LEPEDEA axe directed 
narval to the satellite spin axis. The spin axis of IMP-6 is 
oriented perpendicular to the ecliptic plane to within e few degrees. 
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III. OBSERVATIONS 

With the IMP-6 plasma wave measurements, Gurnett and Shaw 
[1973] have identified a noise band consisting of electromagnetic 
waves trapped in the magnetosphere at frequencies above the local 
electron plasma frequency. This noise band frequently has a very 
sharply defined lower cutoff at the local electron plasma frequency. 

The cutoff in the noise spectrum is caused by the reflection of 
ordinary mode waves at the local electron plasma frequency. An 
example of this noise band is illustrated in Figure 1 which shows 
the spectrum analyzer data for an outbound IMP-6 pass on April 6, 

1972, in the i-ocal-morning sector of the magnetosphere. The 
ordinate for each frequency channel is proportional to the logarithm 
of the electric field strength in the channel. The interval from 
the v ase~ine of one channel to the baseline of the next higher channel 
represents a dynamic range of 100 db. The r ertical bars denote 
the average field strength over an Interval of 81.8 seconds and the 
dots represent the maximum field strength over the same interval. 
Trapped electromagnetic noise of the type discussed by Gurnett and 
Shaw is clearly evident in the 5.62 and 10.0 kHz channels in Figure 1, 
extending throughout the region from the pla sma pause to the mag- 
netopause. The p?asmapause location at about 1710 UT can be identified 
from the rapid decrease in the frequency of the upper hybrid resonance 
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noise band, as indicated in Figure 1, and by the abrupt increase in 
( 

the electric field interference generated by the spacecraft solar 
array in the four lowest frequencies of the spectrum analyzer [see 
discussion by Gurnett and Shaw, 1973] • The magnetopause position 
at 0020 UT is indicated by an abrupt decrease in ^he intensities 
of the trapped electromagnetic noise. Further direct support for 
this identification of the magnetopause position is given later 
with the plasma measurements. The sharp lower cutoff frequency of 
the trapped electromagnetic noise is shown in the frequency-time 
spectrogram of Figure 2. The lower cutoff frequency is first 
clea-ly evident at about 8.5 kHz at 1820 UT and subsequently decreases 
to a minimum value of about 4.0 :Hz at 2005 UT. As shown by Gurnett 
and Shaw [1973J> this cutoff is at the local electron plasma frequency. 
Small fluctuations in the cutoff frequency, presumably due to 
electron density irregularities, are present on a time scale of a 
few minutes. 

Since the cutoff frequency is sharply defined the plasma 
frequency can be determined with excellent accuracy, typically to 
within + 1 percent. The electron density, n £ , can then be car uted 
from the plasma frequency, f p , using the equation (cgs units) 



where e is the electron charge and m g is the electron mass. It 
should be noted that the electron density given by equation (1) is 
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the total electron density, independent of the electron energy, 
provided the electrons are not relativistic. For all magnet ospheric 
regions of interest here relativistic electrons constitute a negligible 
fraction, < 1$, of the total electron density. 

The electron and proton intensities measured with the 
LEPEDEA during this outbound pass are shown in the energy-time 
spectrograms of Plate 1. The ordinate of each spectrogram is 
electron energy in units of electron volts, and the abscissa is 
universal time. The detector responses are color-coded from blue 
to red (low to high detector responses) according to the log scale 
shown on the right side of the spectrogram. The four uppermost 
spectrograms give the proton energy spectra in each of four quadrants. 
From top to bottom the angular fields of view of the detector are 
directed toward the sun, toward local cxawn, away from the sun, and 
toward local dusk, respectively, as indicated by the sketch of the 
viewing quadrant to the left of each spectrogram. A sector of width 
40° and centered in the solar direction has been eliminated from the 
quadrant directed toward the sun. This subsector spans th solar-wind 
ion flow and is processed separately. The fifth panel from the top 
in Plate 1, labeled GM tube, displays the responses of a thin- 
windowed Geiger-Mueller tube which is sensitive to electrone with 
energies E 2 45 keV and protons with energies E ^ 650 keV. The 
responses of this detector axe also sectored into four quadrants and 
are color-coded according to the color of the viewing quadrant 
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sketches shown in the top four spectrograms. The bottom spectrogram 
in Plate 1 gives the electron energy spectrum averaged over a 
complete rotational period of the spacecraft. 

The position of the magnetopause ax about 0020 ITT is clearly 
evident in the observations shown in Plate 1 in the abrupt decrease 
of the GM tube counting rates and the decreases of the average 
electron and proton energies as e spacecraft passes through the 
magnetopause. The bow shock is located at about 0315 UT as indicated 
Ly a rapid decrease in the magnetosheath proton intensities and the 
appearance of solar wind p-"oton fluxes at ~ 1-keV in the sunward- 
viewing quadrant. Proton energy spectra for thin pass are shown 
in Figure 3 at three selected locations between the plasmapause 
(at 1710 UT) and the magnetopause (at 0020 UT). Two distinct 
components are evident in these energy spectra. There is a broad 
energetic component extending from several hundred eV to greater 
than ter keV, with a maximum at several keV. Thib energetic component 
corresponds to that of the quiescent proton ring current [Frank, 

1967a] and 1 b encountered throughout the region between the plasma- 
xuse and magnetopause . Near the plasmapause the spectrum of the 
ring current protons broadens markedly toward lower energies. At 
low energies, less than 100 eV, a second component, rapidly increasing 
in intensity toward lower energies, can be seen in the spectra 
observed at 1822 UT and 2010 UT. This low-energy component is also 
discernible in the pro! or: spectrograms of ilate 1 at about 1720 UT 
through 1900 UT. As shown by the series of spectra in Figure 3 and 
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toy the proton spectrograms of Plate 1 the intensities of this low- 
energy component decrease rapidly with increasing radial distance 
and are completely undetectable after about 2120 UT. 

To determine the contribution of the hot suprathermal ring- 
current protons to the total proton density we have computed the 
densities of protons with energies 88 eV ^ E S 38,000 eV by using 
the direct observations of differential energy spectra of proton 
intensities with the LEPEDSA . The lowe r in tegra tion limit of 88 eV 
used in these computations is chosen to avoid any contribution from 
the aforementioned low-energy component of the proton spectrum. The 
upper integration limit of 38,00 eV is the upper limit of the 
LEFEDEA energy range. Although the energies of the ring current 
protons do extend above 38,000 eV these protons make only a relatively 
minor '':-uribution to the proton densities because of the rapidly 
decreasing intensities with increasing energy. The proton (88 eV < 

ES 38>000 eV) densities have been computed from the spin-averaged 
directional intensities of protons in order to reduce errors caused 
hy anisotropies in the angular distributions. For this pans 
the angle between the spin -xis and the predicted geomagnetic field 
direction varied from about 67° to 75° • The intensity unisotropy can, 
therefore, be determined over a large range of pitch angles . Through- 
out the region of interest, the observed anisotropies over this 
range of pitch angles were less than 25 percent. It Is conservatively 
estimated that the error in determining the proton density, which is 
attributable to this anisotropy, is less than + 5 %. When all sources 
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of error due to intensity anisotropies and detector calibrations are 
considered, it is judged that the determination of proton densities 
with the LEPEDEA is accurate to within + 10 percent. 

The results of the combined plasma wave and LEPEDEA density 
meets u'ements are shown in Figure 4. The solid curve shows the total 
electron density as determined from the plasma frequency cutoff of 
the electromagnetic noise band shown in Figure 2. The dashed curve 
shows the 88 eV < E < 38,000 eV proton density as measured with the 
LEPEDEA* Since the total electron density and the total proton 
density must be equal to maintain charge neutrality, the difference 
between the two number densities in Figure 4 must be attributed to 
protons with energies not in the range 88 eV < E < 38,000 eV. 

Because of the rapidly decreasing proton intensities with high 
energies, this difference cannot be due to protons with energies 
E ^ 38,000 eV. The observed difference between the two densities 
in Figure 4 must, therefore, be due to "thermal" protons with 
energies E < 88 eV. The high-energy tail of these "thermal" protons 
apparently is associated with the low-energy component of the proton 
spectra shown in Figure 3* 

It is evident in Figure 4 that, from 1820 UT to about 2120 UT, 
the low-energy thermal component accounts for a significant fraction, 
20 to 50 percent, of the total proton density. At about 2120 UT, 
the density of the 88 eV S E < 38,000 eV "suprathermal" protons 
increases abruptly. At this same time the low-energy, E < 88 eV, 
component of the proton spectrum disappears below the thresholds of 
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the LEPEDEA measurements (see Figure 3) and the thermal nroton 
density, given by the difference between the two curves in Figure 4, 
goes to essentially zero. In this region beyond 2120 UT, out to 
the limit af our observations at 2250 UT, the proton distribution 
appears to consist almost entirely of suprathermal ring current protons 
with essentially no detectable low-energy thermal component. The 
close agreement between the two density ‘measurements for this and 
several other similar passes can also be taken as a reassuring 
confirmation of the overall accuracy (+ 10^) of the LEPEDEA density 
determination. 
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IV. DISCUSSION 

These simultaneous measurements of the total electron 
densities and the proton (88 eV < F- < eV) densities provide 

direct measurements ui the thermal proton densities in the loir 
density region between the plasmapause and magnetopause. The hi,?h- 
energy tail of the thermal proton distribution can be detected 
also in the lowest energy channel of the LEPEDEA. Since independent 
measurements can be made of the thermal proton density and the 
directional intensity of the high-energy tail of the thermal distri- 
bution, it is possible to make a quantitative determination of the 
temperature of the thermal protons, assuming that these protons have 
a Maxwellian velocity distribution. An example for which we have 
determined the thermal proton temperature using this technique is 
sunmarized in Figure 5 which shows the proton velocity distribution 
function f p (v) computed at two times during the pass analyzed 1 
the previous section. These times, 2100 UT and 2130 UT, are just 
before and just after the abrupt decrease in the thermal plasma 
density at 2120 UT (see Figure 4). The distribution function has 
been computed from the spin-averaged directional intensities and 
with the assumption that the angular distribution is isotropic. 

This assumption is supported by the near isotropy of the observed 
angular distribution at these times. Directional intensity measurements 
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were not obtained at pitch angles near or within the loss cone where 
the actual distribution function is almost certainly not isotropic; 
however, these protons only represent a minor portion of the total 
proton population. At 2100 UT the thermal proton density 
is 0.2 protons cm . Using this density the temperature of the 
Maxwellian velocity distribution has been adjusted so that the high- 
energy tail intercepts the value for the low-energy component of tne 
proton distribution function determined with the LEPEBEA. The best 
fit in this case was obtained with a temperature of T = 80,000 = K. 
Investigation of other passes through the region between the plasma- 
pause and magnetopause indicates that this temperature is not 
unusual. 

Since the ionosphere is thought to be the source of the low- 
energy thermal plasma observed at high altitudes in the magnetosphere 
[Banks, 1972], it is interesting to note that the thermal proton 
temperatures obtained with these IMP-6 observations are much higher 
than typical ionospheric temperatures of about 1000°K to 2000°K. 

Howe ve' our thermal proton temperatures are in close agreement 
with the results of Serbu and Maier l 197^] who reported ion tempera- 
tures of typically 100 # 000°K in the region beyond the plasmapause 
using measurements from the 0G0-5 spacecraft. Although our results 
loifirm hre ion temperatures obtained by Serbu and Maier [1970] our 
ptrotu densities, which are usually about 1 cm or lesB in the region 
' etween the plasmapause and magnetopause, are a factor of 10 to 
100 smaller than those reported by Serbu and Maier, but are in 
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agreement with densities reported by Taylor et al. , [1970] and 
Chappell et al», [ 1971] • 

Abrupt decreases in the thermal proton densities to values 
much less than the supaathermal proton densities, such as that at 
about 2120 ITT in Figure 4, axe another feature frequently observed 
beyon the plasmapause with the IMP-6 instruments. For the example 
shown in Figure 4 this decrease is accompanied by a marked increase 
in the small scale, few minutes, irregularities in the electron 
densities and by an abrupt onset of electrostatic electron cyclotron 
emissions of the type discussed by Kennel et al. , [1970] near half- 
integral multiples of the electron cyclotron frequency. These 
electrostatic electron cyclotron emissions can be seen at about 
1 kHz in Figure 2 commencing at approximately 2127 UT. Similar 
observations during other spacecraft traversals of these regions 
suggest that these electrostatic electron cyclotron emissions are 
quenched whenever the density of low-energy thermal plasma exceeds 
a certain small fraction, on the order of 5 percent, of the supra- 
thermal plasma density. 
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FIGURE CAPTIONS 

Plate 1 Energy-time spectrograms of the electron and proton 

intensities observed during the pass shown in Figure 1. 
The top four spectrograms are the proton intensities 
in four directional quadrants: from top to bottom, look- 

ing toward the sun, toward local dawn, away from the 
sun, and toward local dusk. The fifth panel from the top 
displays the GM tube responses and the bottom spectrogram 
gives the spin-averaged electron intensities. 

NOTE: Plate 1 is to be published in color . 

Figure 1 An outbound IMP-6 pass at local morning on April 6, 1972, 
with trapped electromagnetic noise above the plasma 
frequency (f > f ) in the region between the plasmapause 

£r 

and the magnetopause. Note the abrupt termination 
of the f > f noise band at the magnetopause. 

XT 

Figure 2 Electric field frequency -time spectrogram illustrating 
the sharp lower cutoff frequency of the trapped f > f 
noise band shown in Figure 1. The cutoff frequency is 
at the local electron plasma frequency. The small scale 
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variations in the cutoff frequency ar^ due to irregulari- 
ties in the electron densities. 

Figure 3 The spin-averaged directional proton energy spectra 

measured with the LEPEDEA at three selected points during 
the pass shown in Figures 3 and 2 and Plate 1. The 
suprathermal protons with energies E i 68 eV and peak 
intensities at several keV constitute the quiescent ring 
current . 'The low-energy E < 88 eV component evident in 
the 3.822 UT and 2010 UT spectra is the high-energy tail 
of the thermal proton distribution. 

Figure 4 Comparison of the total electron densities obtained 

from the plasma frequency cutoff and the suprathermal, 

88 eV< E< 38,000 eV. proton densities obtained with 
the LEPEDEA. Note that after about 2120 UT essentially 
the entire proton population consists of suprathermal 
protons . 

Figure 5 The proton velocity distribution functions for Just 

before and Just after the abrupt decrease in the thermal 
proton density at 2120 UT. The thermal proton temperature 
has been determined by fitting the high-energy tail of 
the Maxwellian thermal distribution to the low-energy 
component of the proton distribution determined with the 


LEPEDEA 
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